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,^Zh ■ Abstract 
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We present Mathematical arguments on the existence and uniqueness of weak solutions 
for a special class of non-linear parabolic and hyperbolic equations of the Mathematical- 
■ Physics subject to random initial conditions sampled in L 2 (ft) space (homogeneous statis- 

\0> tical turbulence). In order to achieve such goal, we use the results of compacity on LP 

spaces with weak topologies, instead of the somewhat elaborated results from the methods 
of parabolic variational inequalities often utilized in such mathematical studies. 

o 
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d 1 1 Introduction 

One of the most important problems in the Mathematical-Physics of the non linear diffusion 
and wave damped propagation is to establish the existence and uniqueness of weak solutions 
in some convenient Banach Spaces for the associated non-linear evolution equation (see ref. 
[1], [3], [4], [5]). Another important class of initial- value problems for non linear diffusion 
- damping came from statistical turbulence as modeled by non-linear diffusion or damped 
hyperbolic partial differential equations with random initial conditions associated to Gaussian 
processes sampled in certain Hilbert spaces ([2], [5]) and simulating the Turbulence physical 



phenomena [6]. 

The purpose of this note is to contribute for such mathematical-physicist studies by using 
functional spaces compacity arguments in order to produce proofs for the existence and unique- 
ness of weak solutions for a class of special non-linear diffusion equations on a "smooth" C°° 
domain with compact closure Q C R 3 with Dirichlet boundary condition and initial conditions 
belonging to the space L 2 (£l). This study is presented in Section 2. 

In Section 3, we present similar analysis for to the Wave equation with a non-linear damping 
analogous in its form to the non linear term studied on Section 2 for the diffusion equation. 

In Section 4, we present in some details the solution of the associated problem for random 
initial conditions in terms of our proposed cylindrical functional measures representations pre- 
viously proposed on the literature ([6]) and defined by a cylindrical measure on the Banach 
space L°°((0, T), L 2 (f2)), a new result on the subject. 

Finally in Section 5, we preset a complementary semi-group analyzes on the very important 
problem of anomalous diffusion for random conditions in the path-integral framework. 

2 The Theorem for Parabolic Non-Linear Diffusion 

Let us consider the following non-linear diffusion equation in some strip Q x [0, T] C R A 

= (-AU)(x,t) + A(F(U(x,t)))+f(x,t) (1) 
with the initial and Dirichlet boundary condition 

U(x, 0) = g(x) G L 2 (n) C L 1 ^) (2) 

U(x,t)\ dQ = (3) 

where A denotes a second order self-adjoint uniform elliptic positive differential operator, F(x) 
is a real function continuously differentiable on the extended real line (— oo, +oo) with its 
derivative F\x) strictly positive on (— oo,+oo). The external source f(x,t) is supposed to be 
on the space L°°([0, T] x L 2 (fi)) = L°°([0, T], L 2 (fi)). 

We now state the existence and uniqueness theorem of ours. 
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Theorem 1. In the initial-value non linear diffusion eq.(l) - eq.(3), for any g(x) G L 2 {VL), and 
A = —A (minus Laplacean) there exists an unique solution U (x, t) on L°°([0, T] xL 2 (Q)) satisfy- 
ing this problem in a certain weak sense with a test functional space given by C£°([0, T], H 2 (Q) n 

Proof. The existence proof will be given for a general A as stated below eq.(3). Let {<pi(x)} be 

the spectral eigenfunctions associated to the operator A. Note that each (fi(x) G H 2 (il)f]HQ(il) 

and this set is complete on L 2 (Q) as a result of the Gelfand Generalized Spectral theorem applied 

for A ([1]) since H 2 (Vt) D Hq(Q) is compactly immersed in L 2 {Vt). 

Consider the following system of non-linear ordinary differential equations associated to 

eq.(l)-eq.(3) - the well-known Galerkin system ([1]). 

'r)Tj( n )(r t) \ 

' -,<pAx)) +(AUW(x,t),<p j (x)) I ? ia) 



/ L 2 (n) 

= (V-[(F\U^(x,t)))VU^\x,t)l^))LHn) 

+ (fW(x,t) i <p j (x)) I ? m (4) 

subject to the initial-condition 

n 

UW(x,0) = Y,(9,<PM<i)-<Pi(x) (5) 
i=i 

where the Finite-dimensional Galerkin approximants are given exactly in terms of the spectral 
basis {{pi(x)} as 

n 

U^(x,t) = J2ut l) (t)^(x) (6) 



i=i 

n 



/<">(*,*) = J2(f(x,t),<Pi(x)) LHn) <pi(x) (7) 

i=l 

and ( , )l 2 {q.) is the usual inner product on L 2 (Q). Note that 

n 

^ (B) (^*)| fln = E^ B) (*)W^))| an = 0. (8) 



i=l 

Let us introduce the short notations 



U {n \x,t) = U (n) (9) 



{g,Vi)L*(n) = 9i n) (10) 

(f(x,t), Vl (x)) LHn) = ti n) . (11) 
By multiplying the Galerking system eq.(4) by as usual ([3]) we get the a priori identity 

^ll^ n) lli 2 (n) + (^ {n) ^ (n) )^(n) 

+ / dxF'(U^)(VU^VU (n) ) = (f,U^) L 2 {n) (12) 
Jn 

By a direct application of the Garding-Poincare inequality to the quadratic form associated 
to the operator A, one has that there is a positive constant 'y(Q) such that ([3]) 

(^i/w.i/W)^) > T (n)||i/W||i, (n) . (is) 

This yields the following estimate for any integer positive p 

\j t ( I \U (n) I &<„)) + 7(0) I |t/ (n) 1 1| 2 (n) + 1 1 (^(C/ (n) )) 1/2 VC/^||i 2(n) 

<^{p||/(^*)lli a( n) + V (n) Hi^)}- (14) 

At this point we observe that the square root of the function F'(x) makes sense since it is 
always positive F'(x) > 0. 

By choosing on eq.(14) p big enough and such that a p = 7(fi) — ^ > 0, we have that 

li 2 (n))+«pl|C/ (n) lli 2 (n)<^ll/]|! 2( n)W, (15) 



or by the Gronwall lema, there is a uniform constant M [evenforT = +oo iff e £ 2 ([0, oo], L 2 (f2))] 
such that 

8up(ll^ (B) (*)lli>(n))< ^p {e^ T [/ T rf S ||/||i 2(n) ( S )e + ^+||^||i 2( J)=M. (16) 

0<i<T 0<t<T [ L^O J J 

Note that the Galerkin system of (non-linear) ordinary differential equations for {U- n \t)} 
has unique global solution on the interval [0,T] as consequence of eq.(16) - which means that 
Y.iU^it)) 2 < M - and the Peano-Caratheodory theorem since f(t, •) G L°°([0,T]). 

i=l 

As a consequence of the Banach-Alaoglu theorem applied to the bounded set {U^} on 
L°°([0, T], L 2 (f2)), there is a sub-sequence weak-star convergent to an element (function) U (t, x) e 



L°°([0, T], L 2 (il)) and this subsequence will still be denoted by {U^(x,t)} in the analysis that 
follows. 

A important remark is at this point. Since F(x) is a Lipschtizian function on the closed 
interval where it is defined we have that the set of functions {F(U^} is a bounded set on 
L°°([0,T],L 2 (fi)) if the set {U^} has this property. 

As a consequence of this remark and of a priori inequalite given in details on Appendix A, 
we can apply the famous Aubin-Lion theorem ([1]) to insures the strong convergence of the 
sequence {U^ n \x,t)} to the function U(x,t). As a straightforward consequence of our special 
non-linearity on eq.(l), namely F(x) is a Lipschitizian function, we obtain, thus, the strong 
convergence of the sequence {F(U^ n \x,t)} to the L 2 (f2)-function F(U(x,t)). 

As a consequence of the above made comments and by noting that L 2 (Q) is continuously 
immersed on H~ 2 (Q) and the weak star convergence definition, we have the weak equation 
below in the test space function of v(x, t) G C£°([0, T], H 2 (tt) n H^(Q)) 



rT 

lim / dt 



= lim 

n— s-oo 



(~U {n \ + (AU^,v) mn) - (F(U^),Av) LHn) 

V / L 2 (Q) 

r dt(fw,v) L2{n) (17) 

Jo 



or by passing the weak-star limit 

dv(x, t) 



T 

dt 



o 



U(x,t),—^) +(U(x,t),(Av)(x,t)) LHQ) 



L 2 (Q) 

T 



-(F(U(x,t)),Av(x,t)) LHQ) ] = f dt(f(x,t),v(x,t)) L2{Q) (18) 

J Jo 

This concludes the existence proof or our Theorem 1, since v(0,x) = v(T,x) = 0. 
Let us apply it for a concrete problem of random exponential non-linear diffusion on a cube 
[0, 1/] 3 C -R 3 , mainly for explanation of the above written abstract results 

^jM = V { (| + e -^*>) VU(x, t) } + ^AU(x, t) (19) 

with 

U (x, 0) = g(x) G L 2 {VL) (20) 
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where g(x) are the samples of a stochastic Gaussian process on Q with correlation function 
defining a operator of Trace class on L 2 (Vl). Namely 

E{g(x)g(y)} = K(x, y) E j{L 2 {n)) (22) 

or in terms of the spectral set associated to the Laplacian, i.e., 

g{x)= limf>fW) (23) 

i=0 

with 

E{gig k } = K ik = dxdyK(x,y)(pi(x)ip k (y). (24) 
Jnxn 

As a consequence of the Theorem 1, an explicitly solution of the equation (19) must be 
taken in L°°([0, oo), L 2 (f2)) in the weak sense of eq.(18) and given thus by the Trigonometrical 
series sequence of functions: 

U(x, t, [g]) = n Bm ( £ ' d^(t, [ 9l \) sen f^-x) sen f^-y) sen f^z) } (25) 
"2->oo ^ i,j,k=l ^ ' ^ ' ^ ' J 

where the set of absolutely continuous function {d l 'i ,k {t, [<7o])} on [0, T] satisfy the Galerkin set 
of ordinary differential equations with "random" Gaussian initial conditions 

cf>^(0, [g e ]) = (</(*), sen {^-x^ sen sen (j^ z )) 2 } = »■ ( 26 ) 

Let us now comment on the uniqueness solution problem for the non-linear initial-value 
diffusion eq.(l) - eq.(3). 

In the case under study, the uniqueness result comes from the following technical lemma 

Lemma 1. IfUri) and Up) in L°°([0,T] x L 2 (Q)) are two functions satisfying the weak rela- 
tionship, for any v E C o x> ([0, T],H 2 (Q) n H^tt)) 



dt{ (u {1) - Z7 (2) , -^J ^ + (Z7 (1) - U m , Av) Lm 

- (F(U W ) - F(U {2) ); Av) LHn) } = (27-a) 
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then C/ (1) = U {2 ) a.e in L°°([0,T] x L 2 {Q)) ([5]). 

Let us notice that eq.(27-a) is a weak statement on the assumption of vanishing initial values 
and it is a consequence of 

lim supp ess |y (Px\U {1) (x,t) -Z7 ( 2)(M)| 2 j = 0. (27-b) 

The proof of eq.(27) is a direct consequence of the fact that F(x) is a Lipschitz function 
satisfying an inequalite of the form 

\F(U W ) - F(U {2) )\(x,t) < ( sup F'(x))\U (1) -U (2) \(x,ty, (27-c) 

— 00<2I<+00 

(for a technical proof see the Avner Friedmon book on ref. [5]). However in the case of A = —A 
(minus Laplacian), as stated in our theorem, the idendity eq.(27-a) means that for functions 
of the form v(x,t) = e~ xt v x (x) E C o °°([0, T], H 2 (Q) n H^(Q)) with Av x (x) = -Xv x (x), the 
following identity must holds true 

dt(F(U {1) ) - F(U {2) ),v x ) L 2 m = (28) 

'o 

which means that F(U(i)) = F(C/( 2 )) a.e. and thus, U(i) = C/( 2 ) a.e., since F(x) satisfies the 
lower bound estimative 

\F(x)-F(y)\>( inf \F'(x)\)\x - y\. 

— oo<x<oo 

This result produces a rigorous uniqueness result for A = —A. 

Another important example of non-linear diffusion equation somewhat related to eq.(l) is 
the density equation for a gas in a porous medium with physical saturation. In this case the 
law of conservation of mass 

V-(pV) + ^ = (29) 

— * 

where V is the velocity of gas, p the density of the gas and P is the pressure, together with the 
Darcy's law and the isothermic equation of state, namely 

P = cp 1 (30) 
V = -kVP 



I 

Jo 
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leads to the following non-linear diffusion equation 



dp(x,t) 
dt 



7cW(p 7 Vp)(x,t) 



(31) 



p(x,0) = f(x)eL 2 



(n) 



Since it is showed that for finite volume open, convex with C°°-boundary domain f2 as in 
our study, the function p(x,t) which satisfies eq.(31) should be a bounded function for any 
t > t > (with to fixed) [5], it is clear that one can replace the above written mathematical 
Gas Porous medium equation by a more physical equation of the form taking into account the 
physical phenomena of saturation 



[0,M], where M is the global upper-bound of the function p(x,t). 

One can thus apply our Theorem 1 to obtain an explicitly set of functions {pu\(x,t)} 
converging in the weak star topology of L°°([t , oo), L 2 (Q)) to the solution p^(x,t) of the 
problem eq.(32). It is a reasonable conjecture that in the non-saturation case, one should takes 
the e — > limit on eq.(32), namely: p(x,t) = limp( £ )(x, t); since F( £ )(x) converge to x 7 in 
the C°°-topology of C([0, M]); As a consequence it is expected that a solution for the physical 
equation eq.(31) should be produced. A full technical description of this limiting process will 
appears in a extended mathematical oriented paper. 

3 The Hyperbolic non-linear damping 

We aim on this third section state a theorem analogous to the Theorem 1 of paragraph 2 on 
diffusion but now for the important case of existence of non-linear damping on the Hyperbolic 
initial- value problem on Q x [0, T] with imposed Dirichlet boundary conditions, including the 
new global case of T = +oo (see Section 2) and a damping positive constant v. 




- 7 cV(F (£) (p)Vp)(x,t) 

f(x) e L\n) 



(32) 



with F( e )(x) denoting a differentiable regularizing function of the function x 1 on the interval 



U(x,t) 
~dt 2 



+ (AU)(x,t) 



= —v 



dU(x,t) 
dt 



+ A 
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The L 2 (ft)-initial conditions are given by 

U (x, 0) = g(x) G L 2 (ft) (34-a) 

U t (x,0) = h(x) e L 2 (n) (34-b) 

U(x,t)\ dn = (34-c) 

and now the non- homogenous term f(x,t) is considered to be a function belonging to the 
functional space 

L 2 ([0, T] x ft) n L°°([0, T], L 2 (ft)). (35) 

We have thus the following theorem of existence (without uniqueness) 

Theorem 2. There exists a solution U(x,t) on L°°([0, T] x L 2 )ft)) for eq.(33)-eq.(35) in the 
weak sense with a test functional space as given by C^°([0, T], H 2 {VL) n if^ft)). 

In order to arrive at such theorem, let us consider the analogous of the estimate eq.(12) for 
eq.(33) with = §- t (U {n \x, t)). Namely: 



^ll^ (n) lli^) + v\\u< n) \\h in) 



+ {AUW,UW)v {a) + ii(n£> (n) )) 1/2 w (n) ii| 2(n ) 

= (f,U (n) )mn) (36) 



or equivalently 



ld {\\U™\\h (a) + {AU^\U^) L ^} 

{ll^ (n) lli^) + (^ {n) ^ (n) )L 2 (n)} 

+ \ \(F'(U {n) )) 1/2 VU {n) \\ 2 L 2 (n) < u(AU {n \U (n) ) 



2dt 
+ v 



l "''ll 2 (n) 



(37) 



If one chooses here the integer p such that ^ = v, we have the simple bound 
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As a consequence of eq.(38), there is a constant M such that the uniform bounds holds true 
(even if for the case T = +00 for the case of f G L 2 ([0, 00), L 2 (f2))). 

supess < t < T ||l^||| 2(n) <M (39) 

sup ess <t<T \\U (n) \\l2 (n) < M. (40) 

As a consequence of the bounds eq.(39) - eq.(40), there are two functions U(x, t) and P(x, t) 
such that we have the weak star convergence on L°°([0,T] x L 2 (Vi)) 

weak-star lim (U {n) (x, t)) = U(x, t) (41) 

n— >oo 

weak-star lim (— U (n) {x, t)) = P(x, t). (42) 

n->-oo at 

We have thus that the relationship below hold true for any test function v(x, t) G C£°([0, T] x 
H 2 (Q) fl ifg(f2)) obviously satisfying the relations t>(x,0) = t>(x,T) = At>(x,0) = Av(x,T) = 
v t (x,0) = v t (x,T) = v tt (x,0) = v tt (x,T) = as a consequence of applying the Aubin-Lion 
theorem in a similar way it was used on eq.(17). 

r-T 



r-T 
'0 

with the initial conditions 



Jo 1 \ a 1 / L 2 (n) 

= f dt(f,v) L 2 {n) (43) 
Jo 



Let us show now that 



U (x, 0) = g(x) G L 2 (f2) (44) 

P(i,0)=li(i)Gl 2 (O). (45) 

C7(x,t) = / dsP(x,s). (46) 
Jo 
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Firstly, let us remark that integrating on the interval < t < T the relationship eq.(36), 
one obtains the following estimate 



^(ll^Wll^n)- 11^(0)11^, 

* d8\\Un\\hv)(s)) + \{AU^{t)M n \t)) L ^) 



Jo J 

- ±(Atf<">(0), tf< B >(0)) + f ds\\{F\U^)) 1 ' 2 VU^\ 2 L^ ) 



o 



2 

<^jym\Un) + ^j\\Un(s)\\ 2 ds (47) 
Since the operator A satisfies the Garding-Poincare inequality on L 2 (f2) 

(AU^,U^) LHn) (t) > l(mu {n) \\h {n) (t), (48) 



one can see straightforwardly from eq.(47) by choosing 2p' > \ and the previous bounds eq.(39) 
that there is a positive Constant B such that 

7 ds\\U^(s)\\ 2 L2{n) <B = MT (49) 



/ 

Jo 



which by its turn yields that dUn ^ i s weakly convergent to P(x,t) in L 2 ([0, T], L 2 {Vt)). 

As a consequence of general theorems of Function convergence on space of integrable func- 
tions (Aubin-Lion theorem again ([1])), one has that P(x, t) is the time-derivative of the function 
U(x,t) almost everywhere on f2, since it is expected that U^ n \x,t) should be a strongly conver- 
gent sequence to U{x,t) on the separable and reflexive Banach Space L°°([0,T] x L 2 (VL)) ([1]). 
(See Appendix B for mathematical details). 

4 A Path-Integral Solution for the Parabolic Non-Linear Diffusion 

Let us start the physicist oriented section of our note by writing the abstract scalar Parabolic 
equation eq.(l) - eq.(3) in the integral (weak) form 

U(t) = e~ At U(0) + f ds e-^ A AF(U(s)) (50) 
Jo 

where F(x) is a general non-linear scalar functional such that F: L 2 (Q) — > H 2 (Vt). 
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Let us supose that the initial conditions for eq.(l) are samples of a Gaussian stochastic 
processes belonging to the L 2 (Q) space. For instance, this mathematical fact is always true 
when the processes correlation function defines a integral operator of the trace class on L 2 (Q) 

E{U(0,x)U(0,y)} = K(x,y) (51) 

where 

/ dxK(x,x) < oo (52) 
Jq 

At a formal path-integral method, one aims at compute the initial condition average of 
arbitrary product of the function U(x,t) = U x (t) for arbitrary space-time points. This task is 
achieved by considering the associated characteristic functional for eq.(51), namely 

ZW)\ = ~^-E jexp (i£dt(j(t),U(t)) LHn ^J | . (53) 

In order to write a path- integral representation for eq.(53), we follow our previous studies 
on the subject ([6]) by realizing the random initial conditions U(0) = U average as a Gaussian 
functional integral on L 2 (Q), namely: 

2[/W] = i / dfi[U ]x [ J] dU(x,t) 



xen,ie[o,T] 



U(t) - (e- At U + J^ds e-( f -^AF([/(s))j 
xexp^ dt(j(t),U(t)) L 2 m ^ (54) 



where formally 



= fnW^))) x JjxdyU o (x)K- 1 (x,y)U (y)Y (55) 

By using the well-known Fourier functional integral representation for the delta-functional 
inside the identity eq.(54) ([6]) (see this old idea on Monin A.S. and Yaglom A.M. - "Statistical 
Fluid Mechanics" - Mit Press, Cambridge, 1971, vol. 2)) ([6]), one can re- write eq.(54) in the 
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following form 

z\j{t)\ = { / nw 1 -')) x / dx ^ 

x exp ii jT ^A(t), l/(t) + jf ds e- ( *- s)A AF(f/)(s) 
xexp|-i^ S <A( S )((e-^- 1 e^)A)( S )> L2(n) }. (56) 



L 2 (Q) 



After evaluating the Gaussian cylindrical measure associated to the Lagrange multiplier 
fields X(x,t), one obtains our formal Path-integral representation, however with a weight well- 
defined mathematically as showed in Section 2, as the main conclusion of Sections 2, 3, 4 



Zj { {J ) U L°°((0,T),L 2 (Q)) 



x exp <( — i y (is y rfs' 



^ + AC/-AF(C/))(s); 



L 2 (Q) 



x + AU - AF(U?) (a') 

xexp^y rfs(j(s),f/(s))i 2 (n)^ • (57) 

It is worth re- write the result eq.(57) in the usual Physicist notation of Feynman path 
integrals ([6]) 



S[/0M)] = i{/ D F [U(x,s)} 

Z l U J L JL°°((0,T),L 2 (n)) 

xexp|-iy dsjf rfs^ tfo^^ + 4l7-AF(E/)^(a:,s) 



x ^£ + A U-AF(U^j (x,s')} 



T 



x exp < i / ds I dx j(s,x)U(s,x) > (58) 



5 Random Anomalous Diffusion, a semi-group approach 
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Recently, it has been an important issue on mathematical-physics of diffusion on random 
porous medium, the study of the anomalous diffusion equation on the full space R D but under 
the presence of a positive random potential as modeled by the parabolic quasi-linear equation 
(0 < t < oo) as written below ([7]). 

= -D (-A) a U(x, t) + V(x)U(x, t) + F(U(x, t)) (59-a) 

U(x,0) = f(x) (59-b) 

where U(x,t) is the diffusion field, D the medium diffusibility constant, V(x) denotes the 
stochastic samples (positive functions) associated to a general random field processes with 
realizations on the space of square integrable functions (by the Minlo's theorem - see Appendix 
D), F(x) denotes the problem's non-linearity represented by a Lipschitzian function as on the 
previous sections and, finally, (— A) a represents the effect of the anomalous diffusion of the 
field U (x, t) on the ambient R D where the diffusion takes place and it is represented here by 
a fractional power of the Laplacean operator. The constant a is called here the anomalous 
diffusion expoent. 

In order to give a precise mathematical meaning we will present a different mathematical 
scheme of the previous sections for eq.(59). Let us, thus, proceed for a moment by re- writing it 
formally in the weak-integral form for those initial- values f(x) G L 2 (R D ) as in eq.(54). Namely 

U(t, [V]) = exp[-t(D (-A) a + V)]f 

+ g f <is{exp[-(t - s)(D (-A) a + V)]}F(U(s)) (60) 
Jo 

where we have used a notation emphasising the stochastic variable nature of the diffusion field 
U(x,t, [V]) as a functional of the samples V{x) G L 2 (R D ) associated to our random diffusion 
potential. 

Physical quantities are functionals of the diffusion field U(t, [V]) and after its determination 
one should average over these random potentials samples. The whole averaging information is 
contained in the space-time characteristic functional as pointed out on section 4 

Z[j(x,t)}= [ dfi\V] x expii / dx f dtj(x,t)-U(x,t,[V})\. (61) 
Jl 2 (r d ) {Jr d Jo J 
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Here d/i[V] means the random potentials cylindrical measure associated the V (^-characteristic 
functional 

Z[h{x)]= [ dfi[V}ex P i((h,v) L 2 (RD) ). (62) 
Jl 2 (r d ) 

Note that we do not suppose the Gaussianity of the field statistics on eq.(62). 

In order to write a functional integral representation for the characteristic functional eq.(61) 
as much as similar representation obtained on section 4, let us re-write eq.(61) into the weak- 
integral form as done in eq.(56) 



Z\j(x,t)} 



dfi[V] / D F [X] / D F [U] 

L 2 (R D ) JL 2 (R D ) JL 2 {R d ) 



x exp*{(A, U + I(F(U)) + e- tA f) L2{RD) } 



(63) 



where I(F(U)) and e tA denotes the objects on the right hand-side of eq.(60) and the (formal 
at this point of our study) contractive generator semi-group below 



A = A)(-A)° + V. 



(64) 



At this point one could proceed in a physicist way by re-writing eq.(63) as a path-integral 
associated to a dynamics of three fields, (the well-known Martin-Siggin-Rose component path- 
integral - ref . [6] ) . 



Z\j(x,t)]= f d»[V] [ D F [X] [ D F [U] 
Jl 2 (r d ) Jl 2 (r d ) Jl 2 (r d ) 



! (fl D ) Jl 2 {r d ) 
x expi I [A, U], 



o Hf-A)(-A)«) 
i(-f -A)(-A)«) 

x expig(X,F(U)) L 2 {RD) 



L 2 (R D ) 



expi [X,U], 



V 
V 



(65) 



L 2 (R D ) 



It is worth remark that in the usual case of the cylindrical measures d/i[V] be a purely 
Gaussian measure, formally written in the Physicist notation as 



dfi[V] = D F [V]exp \~ [ dx [ dyV(x)K- 1 (x, 

{ 2 J L 2( R D} Jl 2 {r d ) 



y)V{y) 



(66) 
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with K(x,y) defining an integral operator of the trace class on L 2 (R D ), one can easily make 
a further simplification by integrating out exactly the y(x)-Gaussian functional integral and 
producing the effective quartic interaction term as written below 



L 2 (R D ) 





V 








V 


POD 




/ dt i 


dy i 


Jo J 


R D JO 



L 2 (R D ) 



= exp\~ f dx [ dt f dy ! ds(X(x,t)U(x,t))K(x,y)(X(y,s)U(y,s))\ (67) 

I Z J R D JO JR D JO ) 

Let us call the reader attention that at this point one can straightforwardly implement the 
usual Feynman-Wild-Martin-Siggia-Rosen diagramatics with the free-propagator given explic- 
itly in momentum space by ([6]) 

[d t - Dvi-Arr 1 » • (68) 

tw - D (\k\ 2a ) 

It is important remark that all the above analysis is still somewhat formal at this point of 
our study since it is based strongly on the hypothesis that the operator A as given by eq.(64) 
is a Co-generator of a contractive semi-group on L 2 (R D ) which straightforwardly leads to the 
existence and uniqueness of global solution for eq.(59). Let us give a rigorous proof of ours of 
such self-adjointeness result, which by its turn will provides a strong connection between the 
parameter a of the Laplacean power, related to the anomalous diffusion coeficient, and the 
underline dimension D of the space where the anomalous diffusion is taking place. 

Let us first recall the famous Kato-Rellich theorem on self-adjointness perturbative of self- 
adjoint operators. 

Theorem of Kato-Rellich. Suppose that A is an self-adjoint operator on L 2 (R D ), B is a 
symmetric operator with Dom(5) D Dom(A) and such that (0 < a < 1) 

H-B^IU 2 ^) < a\\A(p\\ L 2 {RD) + b\\(p\\ L 2 {RD) . (69) 

Then, A + B is self-adjoint on Dom(A) and essentially self-adjoint on any core of A. 
In order to apply this powerfull theorem to our general case under study, let us recall 
that the fractional power of the Laplacean operator, (— A) a , has as an operator domain the 
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(Sobolev) space of all square integrable functions ip(x) such its Fourier Transform (p(k) satisfies 
the condition \k\ a (p(k) G L 2 (R D ) and as a core domain the Schwartz space S(R D ) ([8]). 
As a consequence, we have the straightforward estimative for if e S(R D ) 

\\<f\\L2(RD) < || ^ 2a \\ L 2 {RD) \\(k 2a + l)lp\\ L 2 {RD) 

< C(a,D)(\\k 2a <p\\ L i (RD) + ||£|U W ) (70) 

here the finite constant 

C( ^I(TT^ <TO (71) 

if a > ^, a condition relating the "anomalous" diffusion coeficient a and the intrinsic space- 
time dimensionality D as said before on the introduction. Note that for the physical case of 
D = 3, we have that a = § + e with e > 0. 

Let us re-scale the function (p(k) (r > 0) (for j3 > arbitrary) 

<p r (k) = r fi ip(rk) (72) 



We have thus 



\vA\lhrd) = / dkr p <f(rk) 

J R D 

= r^ D M\ LHRD) (73) 



and 



together with 



^||L 2 (^)=r (2/3 - D)/2 ||^|| L2( ^) (74) 



\\K 2 «Cp r \\ L > {RD) = rW- D -^\\K 2 "<p\\ L2{RD) (75) 
Let us substitute eqs.(73)-(75) into eq.(70). We arrive at the estimate 



\<Pr 



LHRD) < C(a,D)\r^- D -^ 2 \\K^D\\ LHRD} 



+ r w- D v 2 \m\ L2{RD) 



(76) 
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or equivalently 



mW{R») < C(a,D)(r%- 2a \\K 2a p\\ LHRD) 

+ r%\\{p\\ L 2 {RD ^j. (77) 



Let us now estimate the function Vip on the space L 2 (R D ^ 



\\ v ¥\\l 2 (r d ) < IMIl 2 cr d )IMIl°°cr d ) 

< II^IIl^r^IMU 1 ^ ) 



< 



V\\ L , {RD) C(a, D) ■ rT- 2 «) | \K 2a (p\ \„ {RD) 
+ (\\V\\ L2{RD) C(a,D) ■ r%)\\<f\\ L 2 {RD) (78) 

Now one just choose r such that 

\\V\\ L2(R n ) xC(a,D)xr( D - 4a V 2 <1 (79) 

in order to obtain the validity of the bound a < 1 on eq.(69) and, thus, concluding that 
(— A) a + V is an essentially self-adjoint operator on S(R D ). So, its closure on L 2 (R D ) produces 
the operator used on the above exposed semi-group Co-construction for V(x) being a positive 
function almost everywhere. This result of ours is a substantial generalization of Theorem X.15 
presented in ref. [8]. Note that the V(x) samples positivity is always the physical case of 
random porosity in Porous medium. In such case, the random potential is of the exponential 
form V(x) = V (exp{ga(x)}); with g a positive small parameter and V is a back-ground 
porosity term. The effective cylindrical measure on the random porosity parameters is written 
as (see eq.(66)) 

d/i[V} D F [a(x)}exp\-l [ dx [ dye ga(x) R-\x, y)e 9a(y) \ (80) 

9«1 [ 2 J R D J R D J 

The above exposed Co-semi group study complements the study on non-linear purely diffu- 
sion made on the previous sections by purely compacity arguments. 

Finally, let us briefly sketch on the important case of wave propagation U (x, t) in a random 
medium described by a small stochastic damping positive L 2 (R D ) function u(x), such that 
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is a Gaussian process sampled on the L 2 (R D ) space. The hyperbolic linear equation governing 
such physical wave propagation is given by (see Section 3) 

?W +Kl )«M = . WM (81) 

U(x,0) = f(x)EL 2 (R D ) (82) 
U t (x,0)=g(x)eL\R D ) (83) 

where 



(84) 



*(^^)=ff(*,»)ejf(L'(tf>)) 

On can see that after the variable change 

U(x,t,[u])=e~^ t -^(x,t,[u]), (85) 

the damped-stochastic wave equation takes the suitable form of a wave in the presence of a 
random potential 

d 2 $(x,t) . A ,„_ . fv*(x)\^, , 

^ ; = -(-Ar$(M) + (^^j $(x,t) 

$(x,0) = /(a;) 

^(x,0) = ^/(x) + (/(x)=^( a ;). (86) 
Equation (86) has an operatorial scalar weak-integral solution for t > of the following form 

S(f) = cos(tVA)f + (87) 
here, the non-positive self-adjoint operator A is given explicitly by 

A = (-AT - (88) 

As a conclusion, one can see that anomalous diffusion can be handled mathematically in 
the framework of our previously technique of path integrals with constraints. 
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APPENDIX A 



In this mathematical appendix, we intend to give a rigorous mathematical proof of the result 
used on Section 2 about Parabolic non-linear diffusion in relation to the weak continuity of the 
non-linearity AF(U) used on the passage of eq.(17) for eq.(18). 

Firstly we consider the physical hypothesis on the non-linearity of the parabolic eq.(l) that 
the Laplacean operator has a cut-off in its spectral range, namely 



AF(U(x,t)) ->• A^F(U(x,t)) 



(A-l) 



where the regularized Laplacean A^ A ) means the bounded operator of norm A, i.e., in terms of 
the spectral theorem of the Laplacean 



XdE s (X) 



(A-2) 



it is given by 



A (A) = / XdE s (X) (A-3) 

J\X\<A 

Let us impose either the well-known Path-integral Sturm-Liouville time boundary conditions 
imposed on the elements of the path- integral domain (eq.(58)), when one is defining the path 
integral by means of fluctuations around classical configurations 



U(x,T) = U(x,0) = 



(A-4) 



It is straightforward to see the validity of the following chain of inequalites with a > 1 



Jo 



dt 



dU n {t) 



dt 



L 2 (Q) 



< Re&l(( AU n (t),U n (t)) - (AU n (0),U n (0))) 



f f <H 
Jo 



A^F(U n (t)) ■ 



dlf„ 



dt 

|2 



<0 + i) + a\ j o \\AF(U n (t))\\ LH n)dt 



+ 





dU n 


2 










dt 





(A-5) 
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In others words 



1 

1 - - 

a 



T 



dt 



dU n 



dt 



<caA 2 - [ \\U n (t)\\ 
Jo 

< CaA 2 TM 



L 2 (a) / 

2 

L 2 (f2) 



(A-6) 



where c = sup ^'(or) | and M is the bound given by eq.(16). 

a<x<b 

Now by a direct application of the well-known Aubin-Lion theorem ([1]), one has that 
the sequence {U n } is a compact set on L 2 (VL). So its converges strongly to U on L 2 (f2) as a 
consequence of the Lipschitzian property of F(x). 

We have either the strong convergence of F{U n ) to F(U), since for each t we have that the 
inequalite below holds true 



dx\F(U n ) - F(U)\ 2 L2(n) (t) < ( sup (F'(x)) 



dx\U n -U\ 2 LHn) (t)) ^0. (A-7) 



SI \— oo<a;<+oo / \J Q 

It is, thus, a imediate result the validity of the weak convergence on L 2 (Vt) used on eq.(18) 



rT rT 

lim / dt(F(U n ),Av) L 2 (n) = / dt(F(U),Av) L 2 {n) . 



(A-8) 
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APPENDIX B 



In this appendix, we show that the function U(x,t) e L°°([0, T], L 2 (f2)) the weak solution 
of eq.(l) in the test space D((0, T), L 2 (Q)) (the Schwart space of L 2 (Q)- valued test functions), 
satisfies the initial condition eq.(2). In order to show such result, let us consider a test function 
possesing the following form for each e > 0, v ( f\x,t) = (fi(x) (x[o,e](t) / e) , where ipi(x) is a 
member of the spectral set associated to the self-adjoint operator A: H 2 (£l) (~]Hq(£1) — > L 2 (f2). 
The notation X[a,b](t) denotes the charachteristic function of the interval [a,b]. 

Since U m (x,t) converges to U(x,t) in the weak-star topology of L°°((0, T), L 2 (Q)) and 
D({0,T),L 2 {Q)) C L\{0,T),L 2 {Q)), we have the relation below for each vf\x,t) holds true 



dt{U m (x,t),vl £ \x,t)) L 2 (n) — )■ / dt(U(x,t),vl £ \x,t)) L 2 (n) (1-B) 



^0 



or equivalently: 



-/ dt(U m (x,t),ipi(x)) L 2 {a) — >- ! dt(U(x,t),ipi(x)) L 2 {n) . (2-B) 
£ Jo £ Jo 



By means of the mean- value theorem applied to both sides of eq. (2- A) and taking the limit 
e — > 0, we have for each i E Z that: 

(U m (x, 0), (pi(x)) L 2 {Q) — > (U (x, 0), (pi(x)) L 2 {n) (3-B) 

as a direct consequence of our choice of the Galerkin elements eq.(6) whose functional form do 
not changes on function of the order m G Z. 

As a consequence of eq.(5) and eq.(6), it yields the result 

(g(x), ^)z, 2 (fi) = {U (x, 0), ipi(x)) L 2 {n} (4-B) 

or by means of Par serval- Fourier theorem 

hm\\U(x,t)-g(x)\\ L 2 (n) = (5-B) 

result, which by its turn, means that the function U(x,t) obtained by means of our compacity 
technique satisfies the initial-condition eq.(2) in the L 2 (f2)-mean sense. 
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APPENDIX C 



Le us show that the functions U(x,t) e L°°([0, T], L 2 (fi)) given by eq.(41) and P(x,t) e 
L°°([0,T], L 2 (Q)) - eq.(42) are coincident as elements of the above written functional space of 
L 2 (Vt) - valued essential bounded functions on (0, T). 

To verify such result, let us call the reader attention that since U m (x, t) is wearly convergent 
to U(x,t) in L°°([0, T], L 2 (f2)), we have that the set {U m (x,t)} is convergent to the function 
U(x,t) as a Schwartz distribution on L 2 {Q) since D([0, T],L 2 (Q)) C L 1 ([0, T], L 2 (fi)). 

This means that 

I7 m (a;,t)-H7(a;,t) in D'([0,T],L 2 (n)). (1-C) 

Analogous result hold true for the time-derivative of the above written equation as a result 
of U (x, t) be a function 

Q-U m (x, t) -U(x, t) in D'([0, T], L 2 {Q)) (2-C) 

By the other side, the set j dUm ^'^ j converges wearly star in L 1 ([0, T], L 2 (fi)) to P(x,t) G 
L°°([0,T],L 2 (fi)) which, by its turn, means that 

^M^POM) in r/([0,T],L 2 (fi)). (3-C) 

By the uniqueness of the limit on the distributional space D'([0, T], L 2 (Q)), we have the 
coincidence of 9U q^ and P(x, t) as elements of D'([0, T], L 2 (f2)). However, P(x, t) is a function, 
so by general theorems on Schwartz distribution theory j 9U q^ j must be a function either since 
L 2 (f2) is a separable Hilbert Sapce. As a consequence we have that dU j£'^ = P(x,t) as elements 
of L°°([0,T], L 2 (Q)), which is the result searched 

dU{x^t) = -p (M) ae in Q 0;T ] xa ( 4 _ C ) 
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Appendix D 

Probability Theory in terms of Functional Integrals 
and the Minlos Theorem — an overview 

In this complementary appendix, we discuss briefly and for the mathematical oriented reader 
the mathematical basis and concepts of functional integrals formulation for stochastic process 
and the important Minlos theorem on Hilbert space support of probabilistic measures. 

The first basic notion of Kolmogorov's probability theory framework is to postulate a con- 
venient topological space (Polish spaces) Q formed by the phenomena random events. The 
choosen topology of Q should posesses a rich set of non-trivial compact subsets. After that, 
we consider the a-algebra generated by all opens sets of this - so called - topological Sample 
space and denoted by J-q. Thirdly, one introduces a regular measure d/i on this set algebra J-q, 
assigning values in [0, 1] to each member of J-q. 

The (abstract) triplet (fi,J-h,rf/i) is called thus a probability space in the Kolmogorov- 
Schwartz Scheme ([10]). 

Let {X a } af z^ be a set of measurables real functions on Q, which may be taken as continuous 
injetive functions of compact support on Q (without loss of generality). 

It is a standard result that one can "immerse" (represent) the abstract space f2 on the 
"concrete" infinite product compact space R°° = Y[ (R-a), where R is a compactified copy of 
the real line. In order to achieve such result we consider the following injection of Q in R°° 
defined by the family {X a } ae A 

I:tt^R°° 

(D-l) 

u ->■ {X a (u)} aeA 

A new a-algebra of events on Q can be induced on Q and affiliated to the family {X a } a&A . 
It is the a-algebra generated by all the "finite-dimensional cylinders" sub-sets, which are 
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explicitly given in its generic formulas by 

^Afin = jw G f2 | Afi n is a finite subset of A 
and [(X a ) (m C J] [a a , 

QgA fln 

The measure restriction of the initial measure /x on this cylinder sets of Q will be still 
denoted by /i on what follows. Now it is a basic theorem of Probability theory that /i induces 
a measure z/°°) on FJ (R) a , which can be identified with the space of all functions of the index 
set A to the real compactified line R; F(A, R) with the topology of pontual convergence. 

At this point, it is straightforward see that the average of any Borelian (mensurable) function 
G(oj) (a Random variable) is given by the following functional integral on the functional space 
F(A,R) 

[ df,(u)G(u) = [ &/<«>(/) ■ G(I-\f)) (D-2) 

It is worth recall that the real support of the measure di/(°°\f) in most pratical cases is not 
the whole space F(A, R) but only a functional sub-space of F(A, R). For instance, in the most 
pratical cases A is the index set of a algebraic vetorial base of a vector space E and FJ (R) a 

aGA 

(without the compactifation) is the space of all sequences with only a finite number of non-zero 
entries as it is necessary to consider in the case of the famous extension theorem of Kolmogorov. 
It turns out that one can consider the support of the probability measure as the set of all linear 
functionals on E, the so called Algebraic Dual of E. This result is a direct consequence of 
considering the set of Random variables given by the family {e* } ae A. For applications one is 
naturally leads to consider the probabilistic object called characteristic functional associated 
to the measure when F(A, R) is identified with the vector space of all algebraic linear 

functionals E alg of a given vector space E as pointed out above, namely 

Z[j]= [ du^(f)eMifU)) (D-3) 
where j are elements of E written as j = Yl x a^ a with {e a } a( zA denoting a given Hammel 

aeA fln 

(vectorial) basis of E. 
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One can show that if there is a sub-space H of E with a norm || \\ H coming from a inner 
product ( , )h such that 

/ du^(f). \\f\\% <oo (D-4) 
J EnH 

one can show that the support of the measure of exactly the Hilbert space (H, ( , ) H ) - the 
famous Minlos theorem. 

Another important theorem (the famous Wienner theorem) is when we have the following 
condition for the family of random variables with the index set A being a real interval of the 
form A = [a, b] 

[ d^°°\f)\X t+h (f) - X t (f)\v < c\h\ 1+£ (D-5) 
Je 

for p, e, c positive constants. We have that it is possible to choose as the support of the 
functional measure the space of the real continuous function on A the well-known 

famous Brownian path space 

/ du,{u)G{u) = [ du^(f(x))G(f(x)) (D-6) 

Jfl JC([a,b],R) 

For rigorous proofs and precise formulations of the above sketched deep mathematical re- 
sults, the interested reader should consult the basic mathematical book on the subject: L. 
Schwartz - "Random Measures on Arbitrary Topological Spaces and Cylindrical Measures" - 
Tata Institute - Oxford University Press, 1973 ([10]). 

Let us finally now give a concrete example of such results consider a generating functional 
Z[J(x)] of the exponential quadratic functional form on L 2 (R N ) 

Z[J(x)]=exp{~ / d N x f d N yj(x)K(x,y)j(y)\ (D-7) 

I ^ Jrn Jrn ) 

and associated to a Gaussian Random field with two-point correlation function given by the 
kernel of eq.(D-7) 

E{v(x)v(y)} = K(x,y). (D-8) 

In the case of the above written kernel be a class-trace operator on L 2 (R N ) (f RN d N xK(x, x) < 
oo), on can represent eq.(B-7) by means of a Gaussian measure with support on the Hilbert 
space L 2 (R N ). This result was called the reader attention for on the bulk of this paper. 
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